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Abstract 

In this paper, we study the geometry associated with Schrodinger operator 
via Hamiltonian and Lagrangian formalism. Making use of a multiplier 
technique, we construct the heat kernel with the coefficient matrices of the 
operator both diagonal and non-diagonal. For applications, we compute 
the heat kernel of a Schrodinger operator with terms of lower order, and 
obtain a globally closed solution to a matrix Riccati equations as a by- 
product. Besides, we finally recover and generalise several classical results 
on some celebrated operators. 
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1 Introduction 

We first introduce a second order differential operator with quadratic potentials 

T = -div(AV) + (Bx, x) + (Cx, V) . 

with A, B and C matrices. From now on, we call the fundamental solution of 
the operator dt + T the heat kernel of T . Let us recall some well-known facts 
for B — 0, when T becomes 

H = -div(AV) + {Cx,V). 

Kolmogorov |K34] considers the following equation 

(d t - d 2 Xl - X!d X2 )u = 0, (x 1 ,x 2 ,t) e R x R x K+ 

to describe the probability density of a system with 2n degree of freedom and 
obtains an explicit fundamental solution by Fourier transform. Hormander [H67 
uses the same method to construct heat kernel for H under a condition imposed 
on A and B which is equivalent to the hypoellipticity of dt + H. Beals |B99] 
sketches a method to find heat kernel for H with A > via a probabilistic 
ansatz. 

The study of the generalised Hermite operator L = — A + (Bx,x) is of 
independent interest. It takes Hermite operator and anti-Hermite operator as 
its typical cases. Hermite operator Lh = —A + \x\ 2 arises from harmonic 
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oscillator and has been studied for quite some time (cf. |EG88) . |Gj87| h while 
anti-Hermite operator Lgh = —A — \x\ 2 arises from anti-harmonic oscillator 
discussed in [RS75] . To the best knowledge of this author, the geometry induced 
by anti-harmonic oscillator was seldom studied. [CFllj makes some effort in 
this direction. They study the geometry of generalised Hermite operator Lgh — 
—A+(Bx, x) with B any real matrix by characterizing the behaviour of geodesies 
when spatial dimension n equals 2. 

Our interests concentrate on the case C = 0, and we study Schrodinger 
operator 

L s = -div(AV) + (Bx, x) (1.1) 

where A is a symmetric positive definite n x n real matrix, B is a n x n real ma- 
trix commutative with A, i.e. AB — BA, div, V, and (•, •) denote respectively 
divergence, gradient and Euclidean inner product. In section 2, we quantita- 
tively study the associated Hamiltonian system for any dimension, from which 
we conclude that the singularities are hyperplanes in phase space. This work 
uniformly generalises the results for B positive definite or low dimensional space 
(cf. |CFli] . [Fllaj or |Fllbj ). Moreover, we formally characterize three impor- 
tant objects, namely geodesic, energy and action function of the Hamiltonian 
system, where geodesic formally means ^-component of the solution for Hamil- 
tonian system. All these quantities are given in closed forms, which will play 
a crucial role in constructing explicit heat kernel of Schrodinger operator in 
section 3. 

A common way to compute the heat kernel of Hermite operator is to use 
eigenfunction expansion (cf. |T93j ). In recent, |CCT06] studies Hamiltonian 
system qualitatively from the view point of conservation law of energy, and 
obtains the heat kernel formulae with B a diagonally positive definite matrix. 
[Fllaj and |Fllbj generalise the results in terms of spectral calculus for B any 
positive semi-dchnite matrix. In section 3, we first use the obtained action 
function and a multiplier technique to construct the heat kernel of Schrodinger 
operator Ls when the coefficient matrices A and B are diagonal. It is worth 
mentioning that the heat kernel has slightly different properties as the normal 
one does, primly because that the Schrodinger operator under consideration is 
not linear in x— variables. For this reason, we address the fundamental solution 
of dt + Ls or dt + L as generalised heat kernel. We close section 3 with the 
computation of heat kernel for A and B non-diagonal case. 

The heat kernel has significance in two areas of applications. In section 4, 
we first apply it to obtain the heat kernel for Schrodinger operator with terms 
of lower order 

L = -div(AV) + (Bx, x) + (/, V) + (g, x) + h (1.2) 

where real matrix A is positive definite and commutative with B, f and g are 
vectors, and h is a real number. The heat kernel of L has an ansatz 

K(x, x°;t) = W(t) exp{(a(t)x, x) + (fi(jt)x, x°) + (j(t)x°, x°) + (fi, x) + (v, x )}, 

(1.3) 

where a, /3, 7 are expected to be symmetric n x n real matrices, fi, v to be 
vectors, and we deduce a system of matrix and scalar differential equations as 
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in [599] 



a = AaAa — 



B + B l 
2 



(1.4) 



$ = 4/3 Aa 
j = /3Ap 

fi = AaAfi — 2a f — g 
v = 2/3 Afi - (3f 

W X W = 2tr(Aa) + (A^fi) - (f,fj,) - h 



(1.5) 
(1.6) 
(1.7) 
(1.8) 
(1.9) 



where the dot denotes 4?, and B l denotes transpose of B. The main difficulty 
is to solve the matrix Riccati equation fj 1 -4[) . which is an equation of funda- 
mental importance in control theory [AFIJ03 . Fortunately, the heat kernel of 
Ls provide us a globally closed solution of matrix Riccati equation (|1.4p . and a 
condition to identify the solution of the scalar differential equation (|1.9I) . Then 
other equations and hence the heat kernel of L can be explicitly computed. Last 
section is devoted to the second areas of applications. We will recover and gener- 
alise several classical results on some celebrated operators, including Laplacian, 
Hermite operator and Ornstein-Uhlenbeck operator on weighted space. 



2 Hamiltonian system associated with Ls 

In this section, we consider Hamiltonian system associated with Schrodinger 



with A and B commutative. 

Geodesies, energy and Hamilton- Jacobi action function are three significant 
objects in Hamilton- Jacobi theory and are of their own interest. We study them 
one by one in the following subsections. 

2.1 Geodesies 

The Hamiltonian function of Ls is defined as its full symbol 



operator 



L s 



div(AV) + {Bx,x) 



and the associated Hamiltonian system is 




(2.1) 



Denoting D := 2A (B + B l 



), one has 



x = — 



2Ai = 2A(B + B l )x = Dx, 
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and D is symmetric following from that A is commutative with B. The geodesies 
x(s) between x° and x in M™ satisfy the boundary value problem 



it = Dx 

x(0) = xq, x(t) = x 



(2.2) 



We start with the case when both A and B are diagonal matrices, and write 
them as follows 



r„2 



A = A a 



B = A = 



-b 2 

U 171+1 



-hi 



(2.3) 



(2.4) 



AB = A a A b 



a 2 b 2 



a 2 b 2 

a m u m 



-a 2 b 2 



(2.5) 



where aj > 0, bj > for j £{!,■■■ ,n} satisfy the condition that: 



Condition (C). For i ^ k and l^i,k^morm + l^i,k^n, afbf ^ a 2 }) 2 ,. 



Putting A\ = diag{ a 2}™ u A% = diag{a 2 }] =m+1 , A\ = diag{^}^ =1) and 



A\ = diag{-6|}™ =m+lj ne has 



A = 



At 



Aa 
71 2 



(2.6) 



B = 



and 



A\ 



(2.7) 



D = AAB = 4 



AtA\ 



A%A\ 



(2.8) 



The solution of the linear system x = Dx is a combination of radical solutions 

We 



{e 2a ibj . } m it {e -2a jbj s } m =ii {cos ^ aj b jS ) }? =m+1 , and {sin(2 0j V)} 
write the coefficients in both block and component forms 



' .j=m+i • 
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Ci = 

c 3 = 



Cn 
C21 

C13 

C23 



c 2 
c 4 



C*12 
p22_ 
Cl4 
C24 



where 



Cn = 



C12 — 



di • • 

Cl, TO +l 



Clr, 



Cln 



c 



13 



Cl,n+1 ' ' ' Cl,n+m 



Cl4 — 



c m,n+l 
Cl,n+m+l 



c m,n+m_ 
Ci 2n 



C21 — 



C22 — 



C*23 — 



Cm+1,1 •' 

C„i 
Cm+l,m+l 



Cm+l,n 



Cm+l,n+l ' ' ' Cm+l,ro+m 



G 



24 



Cn,n+1 
c m+l,n+m+l 



Accordingly, we first write the solution vectors as 



afi (a) = (e 2aibl V- - , e A,m6m »)* ) 

£2(5) = (cos(2a TO+ ib m+ is), • • • , cos(2a„6„s))*, 

x 4 (s) = (sin(2a m+ i6 m+ is), • • • ,sin(2a n 6„s))*, 



then 



Cll C12 C13 C14 
C21 C22 C23 C24 



^i(s) 
z 2 (s) 
^3(s) 
x 4 (s) 



C m +l,2n 
Cn,2n 



C21 



C12 
C22 



C13 
C23 



C14 
C24 



'A\A\ 



A%A\ 



A\A\ 







xi(s) 






x 2 (s) 






X3{s) 


A\A\_ 




_x 4 (s)_ 



C XX A\A\ C l2 A a 2 A\ C W A\A\ C u A a 2 A b 2 
C 2X A\A\ C 22 A a 2 A h 2 C 2Z A\A\ C 2 ±A a 2 A\ 



xi(s) 
x 2 {s) 

xs{s) 
x 4 (s) 
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Dx(s) 



A\A\ 



Cn C12 C13 C14 
C21 C22 C*23 C24 



a^2(s) 

x 4 (s) 



A\A\C XX A\A\C 12 A\A\C 1Z A\A\C x i 
A2A2C21 A2A2C22 A2A2C23 A2A2C24 



x x (s) 

X2{s) 

X3, (s) 



Noting x(s) — Dx(s), and the condition (C) implies 

di = diagjc^}^!, C13 = dia,g{cj, n+ j}jL x , 
C22 = diag{c w }" =m+1 , C 24 = diag{c J>+i }y =m+1 . 

Similarly, that ajb^ > for j E 1, n implies 

C12 = C21 = C14 = C23 = 0. 

Thus, 



(2.9) 



x{s) 



C n C 13 
C22 C24 



xi(s) 

x 3 (s) 
x 4 (s) 



where Cn, C22, C13, C24 are diagonal matrices commutative with A°j and A b j 
for j = 1,2. 

Next, the boundary condition in (|2.2[) will establish C^'s. As before, we 
introduce some notations. 



x i = ^ ) ■ 

( (!) 
= (X\ ',- 



C\\ — (cn, • • • , c mm )*, 

C*22 = (c m +l,m+l, ' ' ' ^ 



.X = [X\ 



m+1 ' 
(1) 

m+1 ' 



T (oh* 

) I 1 

1 X n ^ ) ! 



x 2 = {x\ 

,n+l ; " " " 5 Cm,n+rri) s 
C*24 = (c m +l,m+n+l, " * * j Cn,2nY ■ 



Given a non-singular matrix M, we define -jg :— M 1 iV. 

By boundary condition in (|2.2p . Cn, C22, C13 and C24 satisfy the following 
linear equations 








1 1/7 



cos(2ty/ -A%A%) sin(2t x /— ^2^2). 





C u " 




\4] 




C22 




x 2 




C13 




X X 




_C24_ 




_X2_ 



To move on, we make assumption (*) that sin I 2ty/— A^A\ ) is non-singular 



Indeed, the region for sin fa^—A^A 1 ^) singular consists of countably many 
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hyperplanes in (x, t)— space, thus, it has no contribution to the Hamilton- Jacobi 
action function. In section 3, we will give a formal remark about this point. 







1 in 


,-2ty/A<tA\ 







cos(2i v /^4p|) sm.{2tyJ-A%A\) 



,.() 

•'-2 



x 2 . 



If, 







I n -r 













o 



a;. 



(*), 



in 





I III 





J-n—ni 

































In — m 











I 7) I 












sm(2*V = 3lp|) 

: 
: 

2 

: e - 



X2 - cos(2t v /^4pf)a;^ 



■' 1 



v 7 ^ 



sin(2t v /^4p|) : a; 2 - cos(2ty '- A%A b 2 )xl 

i ™0 



,,o 

•' 2 



e 2 V^j 

/ r ' / if x l 



I'n—m 

We read off the solution 



:(2tv /Z ^p2) a; 2 - cot(2tv / ^Pf)^2. 



C u - 



^22 - I§, 



— X\ ; 

1 e 4t, /^i 



' 1- 



Cl3 

C24 



a U^AlA\ 



(2.10) 



1 



1 



csc I 2*-\/-4!Mo I x 2 - cot ( 2t\/-/L§yl! 



15/15 j x 2 

and we may recover 's in (|2.9[) from C,,, 's in (|2.10j) with 

'Cn,e 



' 2- 



3 

I S~1 n—m 



C j,n+j — yCi3, e 
c i,n+j = ( C24, e!j 



3 

n—m 



j = 1 , m 



(2.11) 



, i = m+l,n 

where e™ denotes a m-dimensional canonical basis vector with j*' 1 component 
one and others zero, and e™r™ is defined in the same way. 



7 



Finally, we conclude the previous deduction as 



Proposition 2.1. Suppose that A, B take the form \2.6\) . j2.7\ ). 

geodesies of Hamiltonian system \2. 1\) that solves boundary problem 
t 7^ J*\ , k 6 N + , j £ m + l,n are given by 



Then the 
with 



x(s) 









C22 



C13 







C24 



xi(s) 
x 4 {s) 



Cnxi(s) 



Ci 3 x 3 (s) 

C 2 iXi{s) 



(2-12) 



where CV,-'s and components therein are identified by (|2.9p - (|2.11| . 



2.2 Energy 



By use of Hamilton- Jacobi theory, the energy is conserved along the geodesies. 
In order to compute such energy and consequent action that both are associated 
with A, we introduce M— inner product (• , • ) M :— (M- , • ) with M a symmetric 
positive definite matrix. Indeed, 



d_ 

ds 



[x(s),x{s)) (Ac 



- (x(s),x(s)) (Aa 
= 2(x(s),x(s)) {Aa) _ 1 - (x(s),x(s)) {Aa y 
= (x(s),x(s)) {Aayl - (x(s),x(s)) (Aa) - 1 
= (Dx(s),x(s)) (Aa) -i - (Dx(s),x(s)) (Aa 
= 0. 



(x(s),x(s)) {Aa 



Const =: 2E 



(x(s),x(s)) {Aayl - (x(s),x(s)) (Aayl 

The main task of this subsection is to find such constant E in terms of boundary 
data. In the following deduction, /(T) denotes spectral calculus of continuous 
function / on the selfadjoint operator T . As we know in the previous subsection, 



x(s) 









C22 



Cl3 







C24 



xi(s) 

X2(s) 
£30) 

x A {s) 



x{s) 



and 



x(s) 



Cn 
C 22 

2y/A a 1 A b 1 



Cl 3 

c 24 



-2^/AjAj 



-2^A°~A\ 





xi(s) 




X2(s) 




X3(s) 




_Xi(s)_ 





~A\A\ 




x\{s) 


'Cn C13 " 




A a 2 A\ 




X2(s) 


C22 C24 




A\A\ 








A%A\_ 




_x 4 (s)_ 
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A direct computation shows 

(i(s),i(a)) (jla) _i 

= A[x 1 {s) t 1 x 2 (s)*, x 3 (s)*, x 4 (s)'] 

<M -CnCis^S 

o -c| 4 4 

-C u Ci Z A\ C? 3 4 

o c 22 c 24 /i^ o 








xi(s) 


C 22 C 24 /1 2 




X2{s) 









— C 22 ^ 2 




x 4 (s) 



and 



(x(s),x(a)) (Aa) -i 
= 4[xi(s)*, x 2 (s)*, x 3 (s 
CftA? 



, X 4 (s)*] 

CiiC\-zA\ 


c 22 c 24 yi 2 ' 



C 22 ^ 2 





C 22 C 24 yl' 


C 24 y1 2 





xi(s) 




a; 2 (s) 




«3(s) 




x 4 (s)_ 



So, 



(±(«),i(a)) (AO) _i - (£(«), a;(s)) (ila) _i 
= 4[xi(s)*, x 2 (s)*, x 3 (s)*, x 4 (s)*] 









-2CuCi 3 ^ 







xi(s) 





— (C 22 + C 24 )A 2 










x 2 (s) 


2C11C13AS 













x 3 (s) 











— (C 22 + C\ 4 )A\_ 




_x i (s)_ 



= -IQx^sfCuC^Alx^s) 

- 4x 2 ( s ) t (C 2 2 + C 2 24 )A b 2 x 2 (s) - 4x 4 ( s ) 4 (C< 2 2 + C 2 24 )A\x 4 ( s ) 
= 4{-4tr(CuC 13 4) - tr[(C| 2 + C| 4 )4]}. 
Hence, 



(x,x). 







2 {-4tr(CiiCi 3 il5) - tr [(C 2 2 + C 2 4 )4] } . 



Making use of CVj's solved previously, we have 
tr (CnC^Al) 

= ( C\i, C13 



A\ 



1 



4 \sinh 2 (2t v /Ip[) 
1 



Xl,Xi 



4 



(2.13) 



„o „,o 



4 \ sinh 2 (2t 1 /Ip[) 



4cosh(2tylp[) ^ y 
2 \ sinh 2 (2tv^lp[) S 
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and 



^2 

y 24 



tr [{C\ 
- ( C24, C24 



Ch) A b 2 ] 



C22 ! C*2 



b ^2-cos(2V5Ip|)£| z 2 - cos(2V=gp|)^ . , t Q . 

yl 2 T^r: / 77. tT s ' / = / ' \ IV 2 X 2i X 2l 



sin(2i v /^4p|) 

-x 2 ,x 2 ) + 



sm(2tyJ-A%A b 2 ) 



yl 2 



An 



\2ty/^AjAj) 



l (2tsJ-A*A\ 



Finally, we conclude the following proposition on the energy 

Proposition 2.2. Suppose that A, B take the form V2.6]) . {2.ty . Then energy of 
Hamiltonian system 12.1]) conforms conservation law along the geodesies \2. 9]) 
with constant E given by 



E = -^{{x,x) {Aa 



2 {-4tr (CiiCia^) - tr [(£. 
A\ 



smh 2 {2t^/AlA\) 



Xi,Xi ) + 



A\ cosh{2t^J A<[A\) 
t sinh 2 (2t v ^i?lfr Xl ' Xl / 



-A b 
71 2 



\2t^AjA T 2 ) 



x 2 , x 2 



C 2 2 2 ) A\] } 
/ A\ 



sinh 1 (2t^J AlA\) 



x 1; x 1 



-A\ 



sm 2 (2t v / -A^A b 2 ) 



-A\ cos(2tV^lp|) ^ y 
sin 2 (2i^Ip|) X2 ' X2 / 



2.3 Action function 

In this subsection, we compute the Hamilton- Jacobi action function S, which 
is a crucial ingredient in the construction of heat kernel. It satisfies Hamilton- 
Jacobi equation (cf. |BGG96a| . |BGG96bj and [BGG97Q 

BS 

— +H(x,VS) = 0. 
ot 

Noting in our case H = — \E, we have S = \ J Edt + c. In the multiplier 
method to be adopted in section 3, the factor i and constant c independent of 
variable t will be absorbed by multiplier and volume element respectively. For 
this reason, we do not differentiate energy from Hamiltonian, and simply define 
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action function as S = — J Edt. Integration by parts shows 

/^cosh [2t^AfA\ 



Ji = 



A\ 



sinh' 



(2t^AjAj 



dt = - 



2 \ A i sinh [2ty/AjAj 



A\ cosh^vMpf 



sinh 2 2ty/A\A\ 



dt 



2 V A ° sinh ht^J A\A\ 



J 3 



-A 



(2ty/-AZA> 



sm 



dt 



_ A b COS \2ty/-A%A% 



-A\ cos {2tyJ-A%A\ 



sin 



dt 




Finally, we have 



Proposition 2.3. Suppose that A, B take the form 112.6]) . Then the 

Hamilton- J acobi action function of Hamiltonian system V2.1\) with boundary 
condition in \2.2)) is given by 




Edt 

^tT cosh (2t^AfA\) 



A i sinh [2ty/AlA\ 

A\ 1 
A i sinh (2ty/AjAj 



cos 



2ty/-A$A* 



J1 2 sm 



(2by/=A% 




b cosh (2t^AlA\\ (j \ 



A t sinh hty/AjAj 

(2ty/-A%A> 



A\ cos 



2 sin {n^-A%A h ^\ 



- A \ 



A 2 sin [2ty/-A%A\ 



X2, x 2 



X 2 , X® 



(2.14) 



3 Heat kernel for L$ 

Given diagonal coefficient matrices, we find heat kernel for L$ via multiplier 
techniques, and discuss its properties similar to the normal heat distribution. 
Heat kernel formulae for non-diagonal coefficient matrices will be given at the 
end of this section. 

3.1 Explicit formulae (diagonal case) 

We start with a basic fact on the action function. 
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Lemma 3.1. Given A, B, energy E and action function S as in \2.6\) , \2. r ty , 
h2.10)) and \2.11\) respectively, the following equalities hold 



(1) \V x S\ 2 A =i(Bx,x)+2E 



(2) tr(AHess(5)) = ^ 



lajbj cosh(2toj6j) 
sinh(2tajbj) 



E 

j=m+l 



2ajbj cos{2ta,jbj 
s'm(2tajbj) 



(3.1) 
(3.2) 



where \-\a '■= y/ (■ , ■ ) a, and Hess(/) denotes Hessian of function f G C 2 

Proof. A direct computation shows that 
(1) 



= 2 



Glf cosh(2ty/A1A\) 
A i sinh(2ty / AlA'l) 
— A\ cos (2t^J-A%A%) 
A 2 sm(2t^J -A%A\) 



X2 



A l sinh(2t^/ 'AIA\) 1 

~EK i x 

^2 sin(2t v /^3TfAf) 



|v*s|i 

= [(^5)', (V^S)'] 













y Ba 5 



= 4 ■ 



A\ cosh 2 : (2iV -4pi) 
4? sinh 2 (2tyip[) 



■Xl,Xi 



/l? cosh(2tyip[) 
^?sinh 2 (2tV^p[) 



Af sinh 2 (2tv / ^P[) 



cos 2 (2<V^Pf) 
^2 sin 2 (2tv / -^^2) 



^2,2:2 



\ ^2 sin 2 (2ty^Ip|) 
= 4( J Bx,ai) + 2£. 
(2) 



~.o o 
x 2 , x 2 



- 2 ■ 



— £l| cos(2^/^gp|) N 
Al S m 2 (2ty/^AjA u 2 ) X2 ' X2 / 



Hess{S) = 2 



Glf cosh(2t- v //l°A;) 
A l sinh(2t v //T°vlf) 



Hence, 



AHess(S') 



/.l Q -l b cosh ( 2 V A ? A ?) 

V 1 1 sinh(2t % /A»A[) 



Thus, 



-A\ cos(2t^ -A%A b 2 ) 
A 2 Bm{2ty/ - A%A\) _ 



/ Z ^cos(2tV=APf) 
v J 2 sin(2t- v /-Ag/i^) 



tr(/4Hess(5)) = ^ 



2ajbj cosh(2tajbj) 
sinh(2ta jbj) 



- E 



2aj6j cos(2tajbj 



sin(2tajbj) 



□ 
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We expect to find the heat kernel of L$ in the following form 

K(x,x°;t) = V(t)e KS( - x ' x °^ 

where the multiplier k is a real number. Making use of (|3.1I) in Lemma 4.1 and 
noticing that 

PK : = (d t +L)K 

= (d t - div(AV) + (Bx,x)) (v(t)e KS{x < X0 ^ 

'V 



= K 



V 
V 



-kE- k 2 \ V x S\ \ - Ktr(AHess(S)) + (Bx, x) 



= K ( y - kE - 4k 2 (Bx, x) - 2n 2 E - Ktr(AHess( 1 5*)) + (Bx, x) 



= 

for t > 0, we choose k 
equation 



— h and let volume element V(t) satisfy transport 

v(ty 



V(t) 



Ktr (,4Hess(S')) . 



(3.3) 



Readers may consult |BGG96a] . [BGG96b and BGG97 for more types of trans- 
port equations. By (|3.2p in Lemma 4.1, we integrate equation (|3. 31) to have 



v(t) = cf[ 



3=1 



1 



sinh (2ta,jbj 



n 

j—m-\-l 



\sin (2tajbj 



(3.4) 



The constant C is determined to normalise the integral in ir— variable of 
the heat kernel K. However, it is easy to see from f|2 . 14[) that the integral is 
divergent if xq ^ 0. K(x,0;t), a propagator from origin to arbitrary point x, 
is called generalised heat kernel as we shall prove next section that it indeed 
has similar properties to the normal heat distribution. We denote K(x,0;t) by 
K(x; t) from now on, then 



K(x;t) = C\[ 



1 



n 

j=m-\-l 



i 



\sin(2tajbj 



(j \sinh(2tajbj 

, 4t ^j = l aj ainh(2ta J b J ) >. X J •> Titf=m+1 aj s in (2t a j 6~ ) >- X j > J 



Making use of / R e x2 dx = y/n, 
/ K(x;t)dx 



( S inh(2ta, 6 7 ) J II 
1 2i6 :) cosh(2taj&j) 



3in(2taj6j) 



n 

3=1 



At Qj smh(2tajbj 



1=1 v °J 



n 

3 = 1 



cosh(2tajbj) 



n 

r n 

n 



1 2t6j cos(2tajbj) 
At a,j shx(2tajbj) 

1 



cos(2tajbj) 
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( h 



By choosing C = n,- =1 — — I , we have arrived at the following proposi- 



tion 



Proposition 3.1. Assume that A and B are diagonal matrices as in 12. 3\) and 
j2.4\ ). Then the heat kernel of the Schrodinger operator L$ = — div(AV) + 
(Bx, x) is 



2tb, V rr ( 2tb 



•Li V a j smh(2ta J 6j ) 7 .J£ +1 V flj sm(2ta j &.,■ ) , 

xe 4t ^j = l a,- sinh(2ta J b J ) - 1 ^^ = 1 aj B i„h(2ta.,' bj ) V x j J ^ 

1 Ap" 2tbj oo a (2fa J 6 J ) (1) 2 2ti,j cos(2taj fcj) ■ (0 ) ■ 3 N 

x e ~4t {L, j=m +1 — a i„(2ta J b J .) t X j ) +2. 3 = m + l — ai „(2ta J b J .) l X j J ) 

if v ™ ^ 1 r (0) r (l) I V" 2t6 J 1 „tO)-(l)N 

X g2t^Z JJ = l .inh(2tajbj) ^ +2 J j = m +l ^ si „ (2t aj b j ) X j X J ^ 

(3.5) 

Remark 3.1. In sake of continuity of a^'s and bj's, heat kernel p.5[) keeps valid if 
the matrix D has multiple eigenvalues or zero eigenvalues. Moreover, condition 
(C) is technical, and can be removed. 

Remark 3.2. Heat kernel p. 51) is complex valued as long as sm(2tajbj) < 0, i.e. 

*e (tg^, f^), fc e N+, j e ^+T^. 

Remark 3.3. Heat kernel (|3.5[) holds if the sub-matrix sin ^2£^/— yijAfj is non- 
singular, which we proposed as an assumption in the previous section. We call 
region fi = {( x, t) <E R" x M+ : t = ^ h ~, k € N+, j € to + l,n} singular 
region and region il c = M™ x M+ \ ft regular region. Briefly speaking, there 
is no geodesic or uncountably many geodesies connecting the given boundary 
points x and x° for t = 2 a\- , while there is a unique geodesic for any given two 
points x and x° ii t =/= 2 ^ b . Here we point out that such singular region has 
no contribution to the Hamilton- Jacobi action function which is regarded as an 
integral of energy in t— variable. 



3.2 Generalised heat kernel 

In this subsection, we show that generalised heat kernel has analogue properties 
to the normal one, that is 

Proposition 3.2. Heat kernel h3. 5}) is said to be generalised in the following 
sense. 

(1) K{x,x°;t) > 0, y(x, x°) e R™ x E n , < t < 1. 

(2) fii i°=0, K(£',t)->1, ast^0 + . 

(3) Fixx° = 0, K(x;t) A S(x), as t -> 0+ 
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where hat denotes Fourier transform on spatial variables, and — > means limita- 
tion in the sense of distribution. 

Proof. (1) It is obvious from formulae (|3.5p for t appropriately small. 
In the rest of this proof, we fix x° = in (|3.5|) and 

2^ \* i r ( 2 Jh 



^)=^nbfe n 



x e 



fJi \a j smh.(2ta j b j )J .J^ +1 \a j sia(2ta j b j )J ^.6) 

iAp m 2tb J °°° h ( 2t °j''j) (- T ,(l)-|2 . 2 ">j c °°( 2t °j''j) ,-(l)^\ 



(2) By properties of Fourier transform 

e^" 2 (0 = e~^ 2 

and 

f(xx)(0 = x- 1 f(x- 1 ^), 

we have 

Kfct) = [ K{x;t)e- 2 ^ x dx 

2tbi \ 2 A ( 2tb 



^ t r § U{ a . sinH L j b j ) ) n 



v aj sinh(2taj6j) / . j-*- \aj sin(2tajbj) 

J — X j — 771 + 1 



f l 2tb -i cosh(2ta 1 -6 1 ) , HK? /,i 

n „ 2^00^(2*0^ (1) 2 



X 



Al l^aj sinh(2ta i 6 j ) / ^ sin(2ta i 6 J -) 



n 



/ 1 2i&,- cosh(2ta j b j )\ 2 _ 4 7r 2 t^- Bi ° h(2f °^ J ' ) 



1 \4-7rf a.,- sinh(2taj6j) 



.7 = 



2tb ■ cosh(2ta^ bj ) '•i 



n 

j— 771+1 



1 2^ cos(2to 3 b 3 ) ^ ^_ 4 . 2 ^^<i^ 2 
4nt ctj sm(2tajbj) 



11 \ v cosh(2ta i & J ) y ' . = ii +1 Vcos(2ta i & J ) 

m 2 a, s„h(2»„ lt ,) 2 " , f aln ( 2to -i,.) 

xTTe b J co S h(ata>;-)S TT g ^ coB(2ta>})«3 

j'=l j=m+l 



as 



f ->■ 0+. 
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(3) We write (HBJ as K(x;t) = K 1 (x;t)K 2 (x;t), where 

^ t ) = f[( sinH 2ta j b j )) A (sin(2taA)) 



• e 
and 



*f 2 (z;t) = - ^-re « . 

(detA) 2 

The proof will be carried out in two steps, 
(i) K 2 {x;t) 4 5(x), as i ^ 0+. 
For any ipeC^°(M. n ), 

K 2 (x;t)ip{x)dx = ip{Q) / K 2 (x;t)dx + K 2 (x;t) [ip(x) - ip(0)] dx. 

(3-7) 

The first term 

/ K 2 (x;t)dx = TT / \=e"^ldxj = ff / e'^Uxj = 1, 

JR" JR flj v 47ri JR 

and the second term 

K 2 (x; t) [<p(x) - <fi(0)} dx = / -i '—e ~ [<p(x) - <f(0)} dx 



k™ (det.4) 2 
/ L (n/I^L^?/) - ^(0)1 -> 0, 



as f ->■ 0+. 

As a result of (13.71) . 



lim / .K^fe f)y(^)rfaJ = v(0) = / ^(a;)w(x)da;, 

i.e. X 2 S, as i -> 0+. 

(ii) K(x;t) A S(x), as i ^ 0+. 

Taking any ip g (R n ), we assume that support of </? is contained in the ball 
{x G 1" : |z| ^ i?}, and that ip is dominated by some constant C everywhere. 
For the sake of 



Ktpdx - tp(0) = I (Kx- l)K 2 ipdx + / K 2 tpdx - tp(0) 
t —> + by checking 



and step (i) lim ( _y + L K 2 (pdx = tp(0), it is sufficient to conclude K — > 5, as 



lim / (Ki - l)K 2 ipdx = 0. 
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Indeed, a variable change y 



makes 



(Ki - i)K 2 <pdx = / {Ki(x;t) - l)K 2 (x;t)<p(x)dx 



A 



(K^V^TrtAiy) - l)ip(V^tAiy)e- 7rM dy 



where 



2ta,ib 



j j 



L \ v sinh(2taj& :) 



n 



j=m+l 



2tajbj 



\sin(2ta,jbj 



-TrTV" 1 2 / It"! »j COeh(2« a ■ bj) N 2 / 2ta j b J .cos(2ta j b J ) \1 

, "L^^^A sinh(2ta ji , 3 ) iJ+L^^l^H S in(2ta 3 b J -) 1 j J 



Noticing that 
obtain 



u cosh(u)— sinh(u) 
sinh(tt) 



0+ and 



ti cos(u)— sin(tt) 
sin(ii) 



— > as t — > + , we 



(Ki - l)K 2 tpdx 



2ta,j bj 



n 



2tcijbj 



-»E"= m+ iVjl 



iin(2ta -b ■) 



1 



e 



as < — > + , which completes the proof. 



dy -> 



□ 



3.3 Explicit formulae (non-diagonal case) 

In order to generalise Proposition 3.1 to non-diagonal case, we hrst rewrite (I3.5[) 
in inner-product form 

/ det tp [2WA a A h ] 

xe -i(( v (2tV^x > x)^ o) _ 1 + ( V (2tV^x^ a ) (AO) _ 1 -2<V'(2 ^ V^^) a: ,x ) (Aa) _ 1 ) 

(3.8) 

where <p(u) = ucoth(u), ^(w) = sin h( M ) • 

Suppose that A and J5 are commutative, then there exists an orthogonal 
matrix P such that PAP 1 = A a and P{B + B^P 1 = 2A b . Putting y = Px, 
system x = 2A(B + B t )x —: Dx becomes y = AA a A h y. Moreover, we have 

detip (2tVA a A b ^ = det if) (tVlPj , 

(<p (W/LM 6 ) V,v) A _ t = ((^ a )"V (2tVlM F ) Px,Px^ 

/ PA _1 PV (2*a/AM&) Px,Px 
^4~V (tv 7 !)) x,x 
= ( tvD ) x, x 
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and similarly, 



ip (WylM 6 ) y°, y ^ ^ = (tp. (ty/D) x°,x< 



(A )" 1 \ V / /A- 1 

According to (|3.8|l . we have arrived one of our main results. 

Theorem 3.1. For given matrix A symmetric positive definite, B a real matrix 
such that A and B are commutative, the heat kernel of Schrodinger operator 

L s = -div(AV) + (Bx, x) 

has the following form 



K(x,x°;t) = (4wty 



det 



detA 



(2.1) 

where D = 2A(B + B% cp(u) = «coth(«), ip{u) = ilTTFR ■ 

Remark 3.4. We use the convention ucoth(tt)| u= o = 1, so Z?coth(_D) = /„ on 
the kernel of D, and so on. 



4 Heat kernel for L 

For the first application of Theorem 3.1, we compute the explicit heat kernel 
of operator L. With mention in introduction, we will adopt the ansatz (jl.3[) 
and solve the associated differential equation system (ll.4l) - (|1.9[) . Of all these 
equations, the most difficult one is the matrix Riccati equation (11.41) . Given the 
results of previous sections, we handle this point in a straightforward way. In 
fact, we have the following 

Theorem 4.1. (Globally closed solution for matrix Riccati equation) 

For matrix A symmetric positive definite, B a real matrix such that A and B 
are commutative, matrix Riccati equation 

B -+- B^ 

a = AaAa (1.4) 

has a globally explicit solution 

a = ^A^ip (tVD) . (4.1) 

Besides, the differential equations system 

(3 = 4[3Aa (1.5) 
7 = (3Af3 (1.6) 
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have explicit solutions 



(4.2) 



where D = 2A(B + B t ), ip(u) = ucoth(u), ^»(u) = 

Proof. The solutions are read out from the Theorem 3.1 since equations (|1.4p ~ 
(|1.6[) are irrelevant to vector /, g and constant /i. □ 

Next, we may integrate equations (|1.7D - (|1.9[) for A and B satisfying condi- 
tion of Theorem 3.1. We point out that for singular matrix _B, solutions are 
formulated in component form. To have concise solutions, we assume that B is 
non-singular and commutative with symmetric positive definite A. 



\i — junction 



1 cosh ( ty/~D 

m = tU"V K- 

2 ^DsmhltVD 



v — function 



D sinh 



(Wd)' 



(4.3) 



(4.4) 



Remark 4.1. Free constants in /i— function and i/— function are absorbed by 
function W(t). 

• W — function 
Making use of 



Aa = 



(AfjL, li ) = 7 \f\ 2 A - 1 - 



coth (ty/D) 



' coth' 



\f\\- 



D 

coth (t\/D 



volume element W(t) satisfies 
-1 



—f,g 

f,9): 



(tVo) 



2(5 + 5*) 



9,9 



W~ X W 



2t 



-tr 



-l/l 



2 

A- 1 



' coth 2 ( t\/D 



2(B + B l ) 



9,9 



h. 



Integration yields 
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where Wo = C 
we have 



det 



sinh (t\fi5) 



and 4>{u) 



Taking f = g = h = 0, 



W = W = V = (47rt) _ i 

Consequently, volume element W is given by 

i 

'detV (ty/Dj 



'det V Uv^D 



det A 



W = (47ri)- = 



det ,4 



(il/li-i+fc)*+W*V5)fl.fl> A * s . (4.5) 



Finally, we yields another main result of this paper. 

Theorem 4.2. For A symmetric positive definite, B non-singular such that A 
and B are commutative, the heat kernel of operator 

L = -div(AV) + (Bx, x) + (/, V) + (g, x)+h 

has the following form 



det 



K{x,x";t) = {4irt)-% 



det A 



-(i\f\ 2 A -l+h)t+(4>(tVD)g,g) A t 3 



kiM A -i-{ — ^- L a,x)+{ 17 --^ 7Wj g,x ) 



x e 



(2.2) 



where D = 2A(B + B*), <p(u) = ucoth(u), ^(u) = <j>{u) 



u— coth(w) 



5 Examples 

For the second application of Theorem 3.1, we demonstrate three examples to 
recover and generalise several classical results on some celebrated operators. 
Notation (p(u) — ucoth(w) and ij){u) = sin ^ M ) keep valid throughout the whole 
section. 

Example 5.1. Generalised Laplacian 

Define generalised Laplacian as Lql = — div(AV) where A is a symmetric pos- 
itive definite matrix. With B = in Theorem 3.1, we have D = 0, ip (t\/~D^ = 

= I n and the heat kernel is given by 

K GL (x,x°;t) = Ml^lM^+M^-V)^) 

(dCt ^ 02 (5.1) 

(47rt)-9 

= Y e 4t 

(det A) 2 
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In particular, taking A = I n , kernel becomes 

x-a; | 2 

K L (x,x°;t) = (Ant)~ 2 e ~ 
which is exactly the Gaussian. 

Example 5.2. Generalised Hermite operator 

Define generalised Hermite operator Lqh = — div(^4V) + (Bx,x) with B > 0. 



A = I n , S = diag{6n« {bj>0) 



By Theorem 3.1, we yield Mchlcr formulae 

L/V™ 2*^00811(2*^) (1) 2 2tbj cosMgtbj) (Q) 2 N (5.2) 

X e 4 * V^-> = 1 sinh(2tt j ) V x j J +2^ = 1 B inh(2t6j ) W'j J J 



- It ( 2 52]=i 



inh(2tt :j ) X j X j 



. yl = diag{a 2 }^ 1 , £? = diag{6 2 }™ =1 (o,- > 0, > 0) 
By Theorem 3.1, the heat kernel for the generalised Hermite operator has the 
form 

-J-/V™ 2tb i I "'' (2 "'i ft 3 ) (-r^h 2 ' V n 2tt, 3 CQ 3 h(2fa 3 b J -) (Q), 2 \ (5.3) 
xe 4*^2^ = 1 sinh(2ta,-b,-) '-J J T^j=l a,; sinh (2t 6^ ) V x j J ^ 

xe «l 2 ^ J= i aj si „h(2ta 3 i, J .) a; 3 ^ ; 

• A > 0, £? = diag{fe 2 }" =1 ( 6j > 0), A and £? are commutative. 
With D — AAB in Theorem 3.1, the heat kernel of the generalised Hermite 
operator is given by 



K(x,x";t) = (47ri)-^ 



det 



det A 



e -i- t (( V (2tV^)x,x) A _ 1 +( V (2tV^)x°.x } A _ 1 -2(4,(2tVAB)x.x°} A _ 1 )_ 

(5.4) 

Example 5.3. Ornstein- Uhlenbeck operator on weighted space 

Define Ornstein-Uhlenbeck operator 

H OU = -div(AV) + Bx ■ V 

with A symmetric positive definite and B any real matrix commutative with A. 
Ornstein-Uhlenbeck 

= -div(AV) + AV(f> • V 
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on Hilbert space L 2 (R™ , e ^dx) is unitarily equivalent to the Schrodinger oper- 
ator 

H = -div(AV) + J|V^|i - idiv(AV0) 
denned on Hilbert space L 2 (R n ,dx): 

= THT- 1 

where T is a multiplication operator defined by 

Tu := e2 u. 

Thus, 

Let (/) take the form (Bx, x) satisfying A\7(f> = Bx, then 
Bx + B l x = V</> = A~ x Bx. 

Hence, 

<f)=^((Bx + S 4 ) a;, ^ = 1 (Bx, x)^ , 

Hess(0) = A _1 .B, 
div(AV^) = tr L4Hcss(0)] = tr(B). 

Consequently, 

H = -div(AV^) + J (B t A~ 1 Bx,x) - ^tr(B). 
By Theorem 4.2 with D = / = g = we have the heat kernel of iJ: 

£T(x,xV) 

/det^(VS*I?)\ 2 

x e -it((v(t^^^)^.4 A - 1 + (v(^^^B^)^^ ) yl - 1 -2(^(tV ^ B7s)x^ ) A _ 1 )_ 

For any 5 e L 2 (R", e"*dx), 

e -«** 5 = Te~ tH T~ l g 

= J K(x,y;t)e~ t ^ 1 g(y)dy 

f 4>(x) 4>{y) 

= / e 2 K(x,y;t)e * g(y)e ^dy. 
Finally, the heat kernel of Ornstein-Uhlenbeck operator on weighted space 
H OU ■ £ 2 (R n ,e~ <B *'''' 4 " 1 ds) -> L 2 (K n , e^'~ : ~^ A ^~dx) 



22 



K u(x,x";t) 



= e 2 K(x,x ;t)e 2 



= (Ant)-? ^ 




x 
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